A Solution of a Tropical Linear Vector Equation* 

Nikolai Krivulin^ 



Abstract 

A linear vector equation is considered defined in terms of idempo- 
tent mathematics. To solve the equation, we apply an approach that 
is based on the analysis of distances between vectors in idempotent 
vector spaces and reduces the solution of the equation to that of a 
tropical optimization problem. Based on the approach, existence and 
uniqueness conditions are established for the solution, and a general 
solution to the equation is given. 

Key- Words: idempotent vector space, linear equation, tropical op- 
timization problem, existence condition. 



1 Introduction 

Many applications of tropical mathematics [H El [3l HI [5l El [7] involve the 
solution of linear equations defined on finite-dimensional semimodules over 
idempotent semifields (idempotent vector spases). One of the equation that 
often arise has the form Ax = d, where A and d are given matrix and 
vector, X is an unknown vector, and multiplication is thought of in terms 
of idempotent algebra. Since the equation can be considered as represent- 
ing linear dependence between vectors, the development of efficient solution 
methods for the equation is of both practical and theoretical interest. 

There are several existing solution methods, including those in [Il|2l[3]. 
In this paper another solution approach is described based on the analy- 
sis of distances between vectors in idempotent vector spaces. The results 
presented are based on implementation and further refinement of solutions 
published in the papers [8l [9l [10] and not available in English. 

We start with a brief overview of preliminary algebraic definitions and 
results. Furthermore, the problem of solving the equation under study re- 
duces to an optimization problem of finding the minimal distance from a 
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vector to a linear span of vectors. We derive a comprehensive solutions to 
the optimization problem under quite general conditions. The obtained re- 
sults are applied to give existence and uniqueness condition as well as to 
offer a general solution of the equation. 

2 Preliminaries 

In this section, we present algebraic definitions, notations, and results based 
on [9l [To] to provide a background for subsequent analysis and solutions. 
Additional details and further results are found in [Tl[2l[3l[ll[5l[6l[7]. 

We consider a set X endowed with addition © and multiplication (g) and 
equipped with the zero and the identity 1. The system (X, 0, 1, ©, (8>) 
is assumed to be a linearly ordered radicable commutative semiring with 
idempotent addition and invertible multiplication (idempotent semifield). 

Idempotency of addition implies that x © x = x for all x G X. For 
any x € X+, where X+ = X \ {0}, there exists an inverse x~'^ such that 

(8> X = 1. Furthermore, the power x^ is defined for any x G X+ and a 
rational q. Specifically, for any integer p > 0, we have x'' = 1, x^ = xP~^x, 
x~P = (x-i)P. 

In what follows, we drop the multiplication sign ® and use the power 
notation only in the above sense. 

The linear order defined on X is assumed to be consistent with a partial 
order that is induced by idempotent addition and involves that x < y if and 
only if X © y = y . Below, the relation symbols and the operator min are 
thought in terms of this linear order. 

As an example of the semifields under consideration, one can consider a 
semifield of real numbers IRmax,+ = (K U {— oo}, — oo, 0, max, +) . 

2.1 Idempotent Vector Space 

Consider the Cartesian power X™ with column vectors as its elements. A 
vector with all components equal to is called the zero vector. A vector is 
regular if it has no no zero components. 

For any two vectors a = (aj) and h = {hi) in X"^, and a scalar x G X, 
addition and scalar multiplication are defined componentwise as 

{a © b}i = a, © 6i, {xa}j = xcj. 

Endowed with these operations, the set X"* forms a semimodule over 
the idempotent semifield X, which is referred to as the idempotent vector 
space. 

Fig. [T] illustrates the operations in the space K^ax + with the Cartesian 
coordinates on the plane. 
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Figure 1: Vector addition (left) and multiplication by scalars (right) in 

''^max,+ • 

Multiplication of a matrix A = (aij) G X"^^" by a vector = (xj) G X" 
is routinely defined to result in a vector with components 

{j4a;}j = diiXi © • • • © ain^n- 

All above operations are componentwise isotone in each argument. 
A matrix is regular if it has no zero rows. 

For any nonzero column vector a; = (xj) G X" , we define a row vector 
x~ = {x~) , where x~ = if 7^ 0, and x~ = <Q otherwise. 

For any two regular vectors x and y, the componentwise inequality 
X < y implies a;~ > y~ . 

If a; is a regular vector, then x x = 1 and xx > I, where / is 
an identity matrix with the elements equal to 1 on the diagonal, and 
elsewhere. 

2.2 Linear Dependence 

Consider a system of vectors ai, . . . , a„ G X™. As usual, a vector b G X™ 
is linearly dependent on the system if it admits representation as a linear 
combination 6 = xiOi © • • • © x„a„ , where xi, . . . Xn G X . 

The set of all linear combinations of ai, . . . , a„ form a linear span de- 
noted by spanjai, . . . , a„} . An example of a linear span in Kmax + is given 
in Fig. m 

A system of vectors ai, . . . , a„ is linearly dependent if at least one its 
vector is linearly dependent on others, and it is linear independent otherwise. 

A system of nonzero vectors ai, . . . , a„ is a minimal generating system 
for a vector b, if this vector is linearly dependent on the system and inde- 
pendent of any of its subsystems. 

Let us verify that if vectors ai, . . . , a„ are a minimal generating system 
for a vector b, then representation of 6 as a linear combination of ai, . . . , a„ 
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Figure 2: A linear span of vectors ai,a2 in K^g^ .^. 
is unique. Suppose there are two linear combinations 

b = Xitti © • ■ • © Xnttn = x[ai © • • • © x'j^ttn, 

where 7^ Xi for some index i = 1, . . . ,n. 

Assuming for definiteness that < Xi, we have b > Xitti > x'-ai. 
Therefore, the term x^aj does not affect b and so may be omitted, which 
contradicts with the minimality of the system ai , . . . , a„ . 

2.3 Distance Function 

For any vector a G X"* , we introduce the support as the index set supp(a) = 
{i\ai 0}. 

The distance between nonzero vectors a,b € X™ with supp(a) = supp(6) 
is defined by a function 

p{a,b)= {b;^ai(Bar%). (1) 

i£supp(a) 

We put p{a, 6) = 00 if supp(a) ^ supp(6) , and p{a, 6) = lifa = 6 = (D. 
Note that in IRmax,+ ' function p coincides for all vectors a, 6 G R"* 
with the Chebyshev metric 

Poo{a,b) = max \bi - ai\. 

l<i<m 

3 Evaluation of Distances 

Let ai, . . . , a„ G X™ be given vectors. Denote by ^ = (ai, . . . , a„) a matrix 
composed of the vectors, and by ^ = spanjai, . . . , a^} their linear span. 

Take a vector d G X"* and consider the problem of computing the dis- 
tance from d to ,4 defined as 

p{A, d) = min p(a, d). 
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Taking into account that every vector a E A can be represented as 
a = Ax for some vector a; G X" , we arrive at the problem of calculating 



p{A,d) = min p{Ax,d). 



(2) 



Suppose d = . Considering that A always contains the zero vector, we 
obviously get p{A, d) = 1. 

Let some of the vectors ai, . . . , a„ be zero. Since zero vectors do not 
affect the linear span A, they can be removed with no change of distances. 
When all vectors are zero and thus A = {0} , we have p{A, d) = if d = (D , 
and p{A, d) = oo otherwise. 

From here on we assume that d 7^ and 7^ for all i = 1, . . . ,n. 

Suppose the vector d = (dj) may have zero components and so be irreg- 
ular. For the matrix A = (aij) , we introduce a matrix A = (aij) as follows. 
We define two sets of indices I = {i\di = 0} and J = {j\aij > (D,i G I}, 
and then deteremine the entries in A according to the conditions 



The matrix A may differ from A only in those columns that have nonzero 
intersections with the rows corresponding to zero components in d. In 
the matrix A, these columns have all entries that are not located at the 
intersections set to zero. The matrix A and the vector d are said to be 
consistent with each other. 

Note that when d is regular, we have A = A. 

Proposition 1. For all x = (xi) G X" it holds that 



Proof. With a regular d the statement becomes trivial and so assume d 7^ 
to have zero components. 

Suppose that p{Ax,d) < 00, which occurs only under the condition 
supp(A£c) = supp(d) . The fulfillment of the condition is equivalent to equal- 
ities anxi © ■ ■ ■ ® ainXn = that must be true whenever = . To provide 
the equalities, we must put Xj = for all indices j such that aij ^ for at 
least one index i with dj = (D. In this case, replacing A with A leaves the 
value of p{Ax, d) < 00 unchanged. 

Since the condition supp(Aa;) 7^ supp(d) imphes supp{Ax) 7^ supp(d) 
and vice versa, the statement is also true when p{Ax, d) = 00. □ 

With the above result, we may now concentrate only on the problems 
when A is consistent with d. 




(D, if i ^ / and j G J, 
aij, otherwise. 



p{Ax, d) = p{Ax, d). 
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In order to describe the solution of problem ([2]), we need the following 
notation. For any consistent matrix A and a vector d, we define a residual 
value 

AA{d) = ^{A{d-A)-)-d 

if A is regular, and A^(<i) = oo otherwise. 

In what follows, we drop subscripts and arguments in A^(<i) and write 
A if no confusion arises. 

Below we find the solution when the vector d is regular and then extend 
this result to irregular vectors. 



3.1 Regular Vector 

Suppose that d is a regular vector. First we verify that the minimum of 
p{Ax,d) over X" in ^ can be found by examining only regular vectors 

Proposition 2. // the vector d is regular, then 

p{A,d) = min p{Ax,d). 

Proof. Take a vector y = Ax such that p{Ax,d) achieves the minimum 
value. If y is irregular and so has zero components, then supp(2/) 7^ 
supp(d) , and thus p{Ax,d) = 00 for all x, including x > d. 

Suppose y = {yi, . . . ,ym)'^ > 0, and assume a corresponding vector x 
to have a zero component, say Xj = . We define the set / = {i\aij > ©} ^ 
and find the number e = mm{aYj^yi\i G /} > 0. 

It remains to note that with Xj = e in place of xj = 0, all components 
of y together with the minimum value of p{Ax, d) remain unchanged. □ 

The next statement reveals the meaning of the residual A = A^(d) in 
terms of distances. 

Lemma 1. // the vector d is regular, then it holds that 

p(A,d) = min p(Ax,d) = A, 

where the minimum is attained at x = A{d~A)~ . 

Proof. Let the matrix A be irregular. Then we have supp(^a;) 7^ supp(<i) 
and p{A, d) = 00. Since, by definition, A = cxo , the statement is true in 
this case. 

Suppose A is regular. Taking into account ([1]), we arrive at an optimiza- 
tion problem to find 

min (d~Ax © (Ax)'d). 

a;ex" 
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Take any vector y = Ax such that a; > 0, and define r = d Ax © 
{Ax)-d > (D. 

From the definition of r, we have two inequalities 

r > d Ax, r > {Ax)~d. 

Right multiphcation of the first inequahty by x^ gives rx^ > d Axx > 
d A. Then we obtain x < r{d^A)^ and (Ax)^ > r^^{A{d- A)')' . 

Substitution into the second inequality results in r > r^^ {A{d^ A)^)^ d = 
r~^A^, and so in r > A. 

It remains to verify that r = A when we take x = A{d^A)^ . Indeed, 
substitution of the x gives r = Ad'^A{d^A)~ © A~^{A{d^ Ay)^d = A. 

Finally note that the above vector x corresponds to the vector y = 
AA{d~A)-eA. □ 

Examples of a linear span A = span(ai,a2) and vectors d in the space 
"^max + given in Fig. [31 




A > 1 




Figure 3: A linear span A and vectors d in K^ax + when A = 1 (left) and 
A > 1 (right). 



3.2 Arbitrary Nonzero Vector 

Now we examine the distance between the linear span A and a vector b 7^ . 
It will suffice to consider only the case when the matrix A is consistent with 
b. 

Theorem 2. Suppose the matrix A is consistent with the vector d 7^ . 
Then it holds that 

p{A,d) = min p{Ax,d) = A, 



where the minimum is attained at x = A{d A) . 
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Proof. Note that for the case when d is regular, the proof is given in 
Lemma[TJ Now we suppose that the vector d ^ has zero components. We 
define sets of indices I = {i\bi = 0} and J = {j\aij > 0,i € /}. 

To provide minimum of p(Ax,d), we have to set xj = for all j J . 
This makes it possible to exclude from consideration all components of d 
and the rows of A with indices in I , as well as all columns of A with indices 
in J. By eliminating these elements, we obtain a new matrix A' and a new 
vector d' . 

Denote the linear span of the columns in A' by A' . Considering that 
the vector d' has no zero components, we apply Lemma [1] to get 

p{A, d) = p{A!, d') = AA'id') = A'. 

Furthermore, we note that the minimum p{A'x',d') is attained if x' = 
A' (d'~ A')~ , where x' is a vector of order less than n. 

The matrix A differs from A' only in that it has extra zero rows and 
columns. Clearly, both matrices appear to be regular or irregular simulta- 
neously. 

Suppose that both matrices are regular. Taking into account that the 
vector d' is obtained from d by removing zero components, we have 

A' = y^{A'{d'-A')-)-d' = ^y{A{d-A)-)-d = A. 

Since the optimal vector x differs from x' only in extra zero components, 
we conclude that p{Ax,d) achieves minimum at x = A{d~A)~ . □ 

In the next sections, we consider applications of the above result to 
analysis of linear dependence and to solution of linear equations. 

4 Linear Dependence 

First we give conditions for a vector d G X™ to be linearly dependent on 
vectors ai,. . . ,an G X'", or equivalently, to admit a representation in the 
form of a linear combination d = xitti © • • • © x„a„ . 

We define the matrix A = (ai, . . . ,a„) and then calculate the residual 

A = y'iAid-Ayyd. 

It follows from Lemma [1] that A > 1. The equality A = 1 means that 
the vector d belongs to the linear span A = spanjai, . . . , a„} , whereas the 
inequality A > 1 implies that d is outside A. In other words, we have the 
following statement. 

Lemma 3. A vector d is linearly dependent on vectors ai, . . . ,a„ if and 
only if A = 1. 
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Now we formulate a criterion that a system ai , . . . , a„ is linearly inde- 
pendent. We denote by Ai = (ai, . . . , aj_i, a,;_|_i, . . . , a„) a matrix obtained 
from A by removing column i, and introduce 

6{A) = min AA,(ai). 

1<J<71 

Lemma 4. T/ie system of vectors ai, . . . , a^, is linearly independent if and 
only if 6{A) > 1 . 

Proof. Clearly, the condition 5{A) > 1 involves that A^^. (aj) > 1 for all 
i = 1 , . . . , n . It follows from Theorem [2] that here none of the vectors 
ai,...,a„ is a linear combination of others, and therefore, the system of 
the vectors is linearly independent. □ 

Let ai, . . . , a„ and bi, . . . , be two systems of nonzero vectors. These 
systems are considered to be equivalent if each vector of one system is a 
linear combination of vectors of the other system. 

Consider a system ai, . . . , that can include linearly dependent vec- 
tors. To construct an equivalent linearly independent system, we imple- 
ment a sequential procedure that examines the vectors one by one to decide 
whether to remove a vector or not. 

At each step i = l,...,n, the vector is removed if Aj (aj) = 1, 

where the matrix Ai is composed of those columns in Ai , that are retained 
after the previous steps. Upon completion of the procedure, we get a new 
system Si, . . . , , where k < n. 

Proposition 3. The system ai,...,afc is a linearly independent system 
that is equivalent to ai, . . . , a„ . 

Proof. According to the way of constructing the system ai,...,afc, each 
vector CLi coincides with a vector of the system ai,...,a„. Since at the 
same time, for each aj , it holds that aj € spanjai, . . . , a^}, both systems 
are equivalent. Finally, due to Lemma HI the system Si, . . . , is linearly 
independent. □ 



5 Linear Equations 

Given a matrix A S '/^"'ixn ^ vector d £ X*^, consider the problem of 
finding an unknown vector a; € X" to satisfy the equation 

Ax = d. (3) 

In what follows, we assume that the matrix A is already put into a form 
where it is consistent with d, and use the notation A = IS.A{d) . 

If a matrix A = (ai, . . . , a„) has a zero column, say , then the solution 
of equation ([3]) reduces to that of an equation that is obtained from ([3]) by 
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removing the component Xi in the vector x together with ehminating the 
column aj in ^4 . Each solution of the reduced equation causes equation ([3]) 
to have a set of solutions, where Xi takes all values in X . 

Suppose that A = 0. In this case, any vector a; G X" is a solution 
provided that b = 0, and there is no solution otherwise. If b = 0, then 
equation ([3]) has a trivial solution x = G , which is unique when the matrix 
A has no zero columns. 

From here on we assume that the vector b and all columns in the matrix 
A are nonzero. 

In the following, we examine conditions for the solution to exist and to 
be unique, and then describe the general solution to the equation. 

5.1 Existence and Uniqueness of Solution 

Application of previous results brings us to a position to arrive at the next 
assertion. 

Theorem 5. If a matrix A is consistent with a vector d 7^ (D, then the 
following statements are valid: 

(1) Equation ([3D has solutions if and only if A = 1. 

(2) // solvable, the equation has a solution 

X = {d-A)-. 

(3) If all columns in A form a minimal system that generates b, then the 
above solution is unique. 

Proof. The existence condition and the form of a solution result from The- 
orem [2j The uniqueness condition follows from representation of a vector as 
a unique linear combination of its minimal set of generators. □ 

We define a pseudo-solution to equation ([3]) as a vector that satisfies the 
equation Ax = AA{d^A)~ . 

Note that the last equation always has a solution that is given by a; = 
A{d^A)^ . Clearly, when A = 1 this pseudo-solution becomes an actual 
solution. Moreover, it follows from Theorem [2] that the pseudo-solution 
provides the minimum distance to the vector d over all vectors in the linear 
span of columns of the matrix A in the sense of the metric p. 

5.2 General Solution 

To describe a general solution to equation we first prove an auxiliary 
result. 
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Lemma 6. Suppose that I is a subset of column indices of a matrix A and 
d G span{aj|i G /} . 

Then any vector xj = (xi) with components Xi = (d^ai)^ if i G I , and 
Xi < {d~ai)~ otherwise, is a solution to equation ([3|). 

Proof. To verify the statement, we first consider that d € span{aj|i G /} C 
spanjai, . . . , a„} and thus equation ^ has a solution xi . For the compo- 
nents Xi of the solution, we can write 

n 

d = Xidi = Xidi © XiCLi. 

1=1 iei" 40/ 

We note that the condition d € span{aj|i G 1} is equivalent to an 
equality 

d = XiCii, 

which is valid when Xi = {d^ai)^ for all i ^ I . 

The remaining components with indices i ^ I must be taken so as to 
satisfy inequalities 

d ^ Xidi ^ Xidi. 

It remains to solve the inequalities to conclude that for all i ^ / , we can 
take any Xi < [d^d-i)^ . □ 

Now suppose that the columns with indices in I form a minimal gen- 
erating system for d. Denote by I a set of all such index sets /. Clearly, 
X 7^ only when equation ^ has at least one solution. 

By applying Lemma [6l it is not difficult to verify that the following 
statement holds. 

Theorem 7. The general solution to equation ^ is a (possible empty) 
family of solutions x = {xi\I G X}, where each solution xj = (xj) is given 
by 

Xi = {d^diy, ifi£l, 
Xi < {d^di)~, ifi^L 

Consider a case when the family reduces to one solution set. Let the 
columns in A are linearly independent. Then there may exist only one 
subset of columns that form a minimal generating system for d. If the 
subset coincides with the set of all columns, then the solution reduces to a 
unique vector x = (d~A)~ . 

Graphical illustration of unique and non-unique solutions to equation 
([3]) are given in Fig. [H 
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Figure 4: A unique (left) and non-unique (middle and right) solutions to 
linear equations in IR^^x + • 
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